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Hysteresis in the Linearized Landau-Lifshitz Equation 


A. Chow^ and K.A. Morris^ 


Abstract —The Landau-Lifshitz equation describes the be¬ 
haviour of magnetization inside a ferromagnetic object. It 
is known that the Landau-Lifshitz equation has an infinite 
number of stable equilibrium points. The existence of multiple 
stable equilibria is closely related to hysteresis. This is a phe¬ 
nomenon that is often characterized by a looping behaviour; 
however, the existence of a loop is not sufficient to identify 
hysteretlc systems, but is defined more precisely as the presence 
of looping as the frequency of the input goes to zero. We 
describe these two approaches to identification of hysteresis 
and demonstrate that both the linear and nonlinear Landau- 
Lifshitz equations exhibit hysteresis. The presence of hysteresis 
in the linear Landau-Lifshitz equation, as well as in a simpler 
system also described here. Indicates that nonlinearity is not 
necessary for hysteresis to exist. 

I. INTRODUCTION 

The Landau-Lifshitz equation describes magnetic be¬ 
haviour within ferromagnetic nanostructures such as 
nanowires [10], [11], [12], [18], [26] and nanoparticles 
[23]. These structures are used in memory storage devices 
such as hard disks or tape recordings. Systems modelled 
by the Landau-Lifshitz equation are generally regarded as 
exhibiting hysteresis [8], [33], [34], [38]. Hysteresis is a 
phenomenon that occurs in many processes [24]. Examples 
include other types of magnetization [6], [16], [24], [29], 
smart materials [30], [31] freezing and thawing processes 
[2], [3], [15], [28], and predator-prey relationships [1], [25, 
Section 1.3]. 

It is difficult to define hysteresis precisely. In much of the 
literature, hysteresis is represented by hysteresis operators 
[1], [7], [8], [21], [32], [35], [36]. The definition of hystere¬ 
sis operator can include operators that are not commonly 
understood as hysteretic [24]. Furthermore, the notion of a 
hysteretic operator is not appropriate for the Landau-Lifshitz 
equation since it requires input-output behaviour to be rate- 
independent; that is, the input-output diagram depends only 
on the values of the input and not on its frequency. 

A common theme in defining hysteresis is that of a 
looping behaviour displayed in the input-output map. The 
shape of the loop can appear quite different under different 
inputs, even for the same system. The Landau-Lifshitz 
equation is an example of this, which is demonstrated 
later in the paper. There are other definitions of hysteresis 
available in the literature, aside from a hysteresis operator. 
One such definition defines hysteretic systems as having 
multiple stable equilibrium points and dynamics that are 
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Fig. 1: Dynamics of a system with two stable equilibria for a 
range of values of the input, and corresponding hysteresis loop. 
Subfigures (a)-(g) show the system variables as the input changes 
with time. The left-hand curve in each subfigure is the energy 
functional with the horizontal axis as the system state and the 
vertical axis as the energy. The right-hand curve in each subfigure 
is the system input-output diagram with input on the horizontal 
axis and output on the vertical axis. 

faster than the rate at which inputs are varied [24]. Another 
definition requires persistent looping behaviour in the input- 
output maps [27]. We describe both approaches and show 
that both the linear and nonlinear Landau-Lifshitz equation 
exhibits hysteresis. 

11. Hysteresis 

Hysteresis can be considered from a dynamical systems 
perspective. 

Definition 1: [24, Definition 3] A system exhibits hys¬ 
teresis if it has 

(a) multiple stable equilibrium points and 

(b) dynamics that are considerably faster than the time scale 
at which inputs are varied. 

The standard definition of stable equilibrium point [37, 
Definition 3.2] with zero input is used here. Condition (b) is 
relative to the speed at which a controlled input is changed. 
In many cases, hysteresis is present but is rate-dependent 
[24]. 

Definition [T] is related to the looping behaviour often 
associated with hysteresis. Consider a system with two 
















stable equilibria and suppose the system is initially at the left 
equilibrium (Figure [T^. If the input is increased, the system 
will tend to stay in equilibrium (Figure with only a small 
move upward along the hysteresis curve. When the input 
is increased enough such that the equilibrium disappears, 
the system moves to the right equilibrium (Figure [^. This 
corresponds to moving along the steepest portion of the 
hysteresis loop. The system stays at the right equilibrium 
(Figure [Tdj i. If the input decreases enough so that the 
right equilibrium disappears (Figure [T^, the system moves 
back to the left equilibrium (Figures[lf| [T^ . For systems 
which move to equilibrium faster than changes in the input, 
the transition from one equilibrium to the other is nearly 
instantaneous and the system behaviour can appear rate- 
independent. For further discussion, see [5], [24]. 

However, the existence of a loop in the input-output map 
is not sufficient to define a system as hysteretic. Consider 
the simple example of a damped second-order system 

y{t) + cy{t) + ky{t)=u{t). (1) 

Writing Q in first-order form with u(t) = 0 and [xi, ^ 2 ] = 
\y, y\ leads to 

xi{t)=X 2 (t) (2a) 

X 2 {t)= —cx 2 {t) — kxi{t) (2b) 


Analyzing the unforced system, the only equilibrium point 
is (0,0) and the eigenvalues are 


Ai,2 — 


—c ± — 4fc 


( 3 ) 


If both eigenvalues have nonpositive real part, then the 
equilibrium point (0,0) is stable. 

For arbitrary initial conditions, j/(0) = j/o and y(0) = j/i, 
the solution to Q with u{t) = sin(a;f) is 

(-Ai - A2a;^) yp + uj 
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where Ai, A 2 are defined in ([^. It follows that 
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If both the real parts of Ai and A 2 are negative, then y{t) — > 
0 as w —0, / — 00 , which means regardless of the initial 
condition, the steady state is zero as w —>■ 0 and f —^ 00 . 

Since Q has exactly one equilibrium, the system does 
not exhibit hysteresis according to Definition [ijFigure 
depicts the input-output curves of Q, with c = 15, /c = 1, 
initial condition (?/(0), y(0)) = (0,0) and with input u{t) = 
sin(a;f) for different frequencies ui. For large oj there is a 




(a) a; = 1 
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(c) uj = 0.01 (d) u} = 0.001 

Fig. 2: Input-output curves for equation (111 with c = 15 and 
k = 1. The initial position is {y, y) = (070) and the input is 
u(t) = sin(a;/:) for various frequencies uj. The uncontrolled system 
has one equilibrium point. As a; —>■ 0, the loops disappear. This 
system does not display hysteresis. 


loop in the input-output diagram. However, as uj approaches 
zero, the loop vanishes. This illustrates a second definition 
of hysteresis. 

Definition 2: [27] A system exhibits hysteresis if a non¬ 
trivial closed curve in the input-output map persists for 
a periodic input as the frequency component of the input 
signal approaches zero. 

For more details on Definition see also [20]. Note that a 
closed curve is said to be non-trivial if it has a non-empty 
interior. For example, a constant function is a trivial closed 
curve. 

Definitions [T] and are complementary, not contradictory. 
Definition [T] refers to the internal behaviour of the system 
state while Definition refers to the system input-output 
behaviour. Note that equation ([T]) satisfies neither definition. 

Modify model Q to obtain the nonlinear system 

y{t) + cy{t) + k {y{t) - y^{t)) = u{t). (4) 

Setting c = 15 and k = —1, the input-output curves 
with zero initial condition and periodic input are shown in 
Figure It is clear from the figure that even for small ui a 
loop persists, which by Definition indicates the presence 
of hysteresis. Observe also the differently shaped hysteresis 
loops for different inputs, which illustrate the complex 
behaviour of hysteresis, even for this second-order system. 

System (j^ also satisfies Definition [T] Rewriting (|^ with 
no input u. 


Xl{t) = X2{t) 

X2{t) = —CX 2 {t) — kxi{t) + kx\{t). 



















The equilibrium points are (0, 0) and (±1,0). The eigenval¬ 
ues of the system linearized about (±1, 0) have negative real 
part and hence (±1,0) are stable equilibrium points. This 
satisfies the first condition in Definition [T] (The equilibrium 
(0, 0) is unstable since one of its corresponding eigenvalues 
has positive real part.) Figure demonstrates that the 
dynamics in equation 0 are independent of the rate at 
which inputs are varied given the range of frequencies. This 
satisfies the second condition in Definition [T] A similar ex¬ 
ample for the displacement of a magnetic beam is presented 
in [24, Example 4]. 



(a) oj = 1 
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Fig. 3; Input-output curves for the nonlinear system in 0 with 
c = 15 and k = —1. The initial position is {y,y) = (0,0) and 
the input is u{t) = sin(a;i:) for various w. 



Fig. 4; Input-output diagrams for the nonlinear system in 0 with 
c = 15, A: = —1. The initial position is {y,y) = (0,0) and the 
input is u{t) = sin(a;f) for cj < 0.0001. At these low frequencies, 
the system exhibits rate independence. 


Consider again the linear equation Q but with fc = 0; 
that is 


and its corresponding unforced system in first-order form is 


Xi{t) = X 2 {t) (6a) 

X2{t) = -CX2{t). (6b) 


The unforced system (j^ has an infinite number of equilibria, 
(a, 0) where a is any constant. The eigenvalues of 0 are 
0 and —c and so the equilibrium points are stable if c > 0. 

For arbitrary initial conditions (j/(0),i)(0)) = (2/o,yi), 
the solution to 0 with u(t) = sin(a;f) is 


y { t ) = 2/0 ± ^ 
± 


^ ^ ± ± 
— cos(wf) — 
wc (w^ ± c^) 


It follows that 


lim y{t) = yo ± — (1 - e "=*) 

^->■0 c 

and hence 

lim y{t) = yo + —. 

t^oo c 
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The steady-state of 0 as w —i 0 depends on the initial 
conditions. Input-output curves are displayed in Figure 
with c = 15, zero initial condition and u{t) = sin(u;f) for 
various values of the frequency uj. The figures indicate that 
a loop persists as uj approaches 0. 

Note that the loop is a smooth circle, and does not 
display the same rapid transition as in Figures 3d and 4. The 
continuum of equilibrium points, rather than disconnected 
equilibrium points, may be the reason for the absence of 
sharp jumps. We will see in the next section that hysteresis 
in the Fandau-Fifshitz equation behaves similarly. 
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Fig. 5: Input-output curves for the linear system in 0. There 
are an infinite number of equilibria, (a, 0). The initial condition is 
(y, y) = (0,0) and the input is u{t) = sm{ujt) for various uj. 


y{t)+cy{t) = u{t). 


( 5 ) 


































III. LANDAU-LIFSHITZ EQUATION 

The Landau-Lifshitz equation is a coupled set of three 
nonlinear partial differential equations, 

rh = m X m.a:x — x (m x m^x) (7a) 
ma;(0, t) = mx{L, t) = 0. (7b) 

m(a;,0) = mo(a;) (7c) 

where x indicates the vector cross product, L > 0 is length 
and z/ > 0 a material parameter. The solution to ([^ is 

m(a;,f) = (mi (x, f), 7712 ( 0 ;, <), 7713 ( 2 ;, Q) 

for a; G [ 0 , L]and 

m X m.xx = ( m 2 'm!^ — 

—mimg + mim'^ — m2m"). 

The prime notation indicates differentiation with respect to 
X. Solutions to (j^ are defined on = £ 2 ([ 0 ,T];K^). 
Existence and uniqueness of solutions to 0 is shown in 
[9]. Equation 0 has been investigated in [4], [9], [19, 
Section 6.3.1]. 

On a molecular level, ferromagnets are divided into 
regions called magnetic domains. Each domain in a fer- 
romagnet is magnetized to the same saturation. Mg. Mathe¬ 
matically, this means ||m(a;,t )||2 = Mg where || • II 2 is the 
Euclidean norm. In much of the literature. Mg is set to 1 
and the same convention is followed here; that is. 


|m(o;,f)||2 = 1 


(7d) 


and the initial condition ||mo (a ;)||2 = 1 - The nonnegative 
parameter 1 / in ( |7a) i is the damping parameter. It depends 
on the magnetization saturation and hence on the type of 
feiTomagnet. 

Theorem 3: The set of equilibrium solutions to the 
Landau-Lifshitz equation ( [7a| ) with Neumann boundary con¬ 
ditions (7bI is m(a;) = a where a G ||a ||2 = 1. 


Proof: It is shown in [19, Theorem 6.1.1] that the equilib¬ 
rium solutions to 01 ) are 


m(a;) = acos(fca;) -)-bsin(fca;) 


where /c G K is a constant, a = ( 01 , 02 , 03 ) and b = 
( 61 , 62 ,^ 3 ) for ai,bi G K, * = 1,2,3 constants with 
||a ||2 = 1, ||b ||2 = 1 and a"^b = 0. Applying the boundary 
conditions in ( [7b] i implies k = 0. □ 

Note that, as for the second-order linear example in 
0 , the Landau-Lifshitz equation has an infinite number of 
equilibria. The set of equilibrium points of 0 is 

E ={a = (oi, 02 , 03 ) G : 

oi, 02, 03 constants and a"'"a = 1}. (8) 


Theorem 4: [19, Proposition 6.2.1] Each point in the 
equilibrium set 0 is stable in the E^'Horm. 

If we introduce the notion of an equilibrium set, then E 
can be shown to be asymptotically stable [13]. 


Since the Landau-Lifshitz equation has an infinite number 
of stable equilibria, Definitionl^suggests it exhibits hystere¬ 
sis. 

Consider the Landau-Lifshitz equation with input u(f): 

rh = m X nixx — x (m x rrixx) + u(f)- 

The input is chosen to be u{t) = (0.001 cos(tLif), 0,0); that 
is, there is a time-varying magnetic field in the mi direction 
and no magnetic field in the 777.2 and 777 . 3 . 
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Fig. 6: Hysteresis loops for mi(x,t) of the (nonlinear) Landau- 
Lifshitz equation with x = 0.6 and u = 0.02. The input is u{t) = 
(0.001 cos(cL;f), 0, 0). The initial condition is mo(a:^) = (1,0,0). 


The hysteresis loops for mi(Xjt) with x fixed is illus¬ 
trated in Figure]^ The initial condition is mo(^) = ( 1 , 0 , 0 ), 
77 = 0.02 and L = 1. It is clear from Eigure|^that the input- 
output curves exhibit persistent loops as the frequency of the 
input approaches zero. In a similar manner, m 2 and 7713 also 
exhibit persistent loops. Thus, these simulations indicate that 
the Landau-Lifshitz equation satisfies Definition 

Note that the hysteresis loops in Eigure have the same 
appearance as those in Eigure for the second-order linear 
example (0. In both examples there is a continuum of 
equilibrium points. 

To obtain the linear Landau-Lifshitz equation, equation 
( f7a| ) is rewritten in semilinear form, 

— = vmxx + mx nixx + i^\\mx\\lni, (9) 

using equation ( [7d| and properties of cross products. Then 
substitute m{x,f) = a + z{x,t) where a G £' is an 
equilibrium of 0 and z G £2 ^ small perturbation into 

0 - Defining 


Az = lyZxx + a X Zxx 













with domain 


D{A) = {z:zeCl, z^e z^^ G Cl, 

Z^{0,t) = 0 = Za:{L,t)}, 

the Landau-Lifshitz equation linearized about an equilib¬ 
rium a is 


z = Az, z(0) = Zq (10a) 

z^{0, t) = 0 = z^{L, t). (10b) 
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The operator A can be shown to generate an analytic 
semigroup [13]. Therefore, (lOi is a well-posed system. 

Theorem 5: Any constant c G is a stable equilibrium 
of ([^. 

Proof: A brief outline is given here; details can be found 
in [13]. Since A generates an analytic semigroup, the 
spectrum determined growth assumption is satisbed and so 
the eigenvalues of A determine the stability of the linear 
system ( [T0| ) [17, Section 5.1], [22, Section 3.2]. 

It is clear that any constant function c is an equilibrium 
of (10 1 . Let A G C. The eigenvalue problem of is 
Av = Av and boundary conditions V2:(0) = Vx{L) — 0 . 
where v G £ 2 - Solving, the eigenvalues of (10 1 are Ai = 0 
and for n G Z, 


-1-.- 

-(1 -f 2n)'^Tr‘^v 

, ,(l + 2n)27r2 

‘2 

£2 

£2 

As = 

— (1 -1- 2n)^7r^z/ 
£2 

1 



(2n)^7r^ 

'4 

£^ 

£2 ’ 

As = • 

— (2n)^7r^i/ 

£2 



Since all the eigenvalues have nonpositive real part, the 
equilibria of ( [T0| ) are stable. □ 

Based on the result in Theorem Dehnition [T] indicates 
that the linear Landau-Lifshitz equation exhibits hysteresis. 
We now look at periodic inputs to determine whether 
persistent loops exists and Dehnition is satisbed. Again, 
ly = 0.02, L = 1, and the same periodic input is applied 
as for the nonlinear Landau-Lifshitz equation. Figure [7] 
shows the input-output curves for the linear Landau-Lifshitz 
equation with a = (1,0,0) and initial condition zq{x) = 
(1,0,0). A loop does persist as the frequency of the input 
approaches zero. Similar plots are obtained when control is 
on the second and third components. The hysteresis loops 
in Figure [7] are similar in shape to the nonlinear Landau- 
Lifshitz equation depicted in Figure and the example 
in each of which have a continuum of equilibria. 


IV. CONCLUSIONS 

The results here indicate that existence of multiple stable 
equilibria is crucial for systems to be hysteretic. The linear 
and nonlinear Landau-Lifshitz equations exhibit hysteretic 
behaviour. Hysteresis was also shown in a second-order 
linear differential equation. Since hysteretic behaviour was 


15 

10 

5 



- 5 . 

- 10 . 

-15 1 . . . . . 

-1 -0.5 0 0.5 1 

“1 X1( 

(c) UJ = 0.01 



(d) UJ = 0.001 


Fig. 7; Hysteresis loops for zi{x, t) of the linear Landau-Lifshitz 
equation with x = 0.6 and v — 0.02. The linearization is at 
a = (1,0,0). The input is u(f) = (0.001 cos(u;f),0,0) and the 
initial condition is zo( 2 ;) = (1,0,0). 


displayed by several linear systems, this suggests that mul¬ 
tiple stable equilibria, not nonlinearity, are the key factors 
leading to the display of hysteresis in systems. As well, 
an absence of jumps in the hysteresis loops were observed. 
Further investigation is required to determine whether this 
is always the case for systems with a zero eigenvalue. 

In [14], the authors discuss controlling hysteresis in 
the Landau-Lifshitz equation; that is controlling between 
equilibrium points. 
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